Localization in momentum space of ultracold atoms in incommensurate lattices 
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We characterize the disorder induced localization in momentum space for ultracold atoms in one- 
dimensional incommensurate lattices, according to the dual Aubry- Andre model. For low disorder 
the system is localized in momentum space, and the momentum distribution exhibits time-periodic 
oscillations of the relative intensity of its components. The behavior of these oscillations is explained 
by means of a simple three-mode approximation. We predict their frequency and visibility by using 
typical parameters of feasible experiments. Above the transition the system diffuses in momentum 
space, and the oscillations vanish when averaged over different realizations, offering a clear signature 
of the transition. 



PACS numbers: 03.75.Lm, 03.75.Kk 
I. INTRODUCTION 



Disorder is among the most intriguing and ubiquitous 
aspects of condensed matter pQ . The prediction of local- 
ization induced by a random disorder in a periodic lat- 
tice dates back to the seminal work by Anderson [2] . Re- 
cently, ultracold atomic gases have demonstrated to be an 
extremely versatile tool to explore the effects of disorder, 
owing to the great tunability of the system parameters 
and geometrical configurations [3H7]. Notably, Anderson 
localization of matter waves has been observed both for 
correlated (speckle) disorder [4 and quasi-periodic op- 
tical lattices [5 , the latter case realizing the so-called 
Aubry- Andre ( AA) model [8j [9] . Localization effects in 
the quantum dynamics of one-dimensional lattice models 
have attracted a large interest also in the recent theo- 
retical literature [T0HT2] . The AA model has the pecu- 
liar property of exhibiting a transition, already in one 
dimension and both in real and momentum space (dual- 
ity), from extended to exponentially localized states as 
the disorder strength is increased above a critical value 
[9] [13] [14]. So far the evolution of wave packets in the A A 
model has been investigated mainly in real space, looking 
for signatures of the transition from ballistic spreading to 
sub-diffusion and localization, both in theory p~0j[l5] and 
experiments [5j [14] . 



Here we focus on the dynamics of the momentum dis- 
tribution and identify measurable effects of the transition 
from diffusion to localization in momentum space. This 
is relevant for current experiments with ultracold atoms, 
where the momentum distribution is accessible via time 
of flight measurements and, typically, with an higher ac- 
curacy than in real space. In addition, it provides com- 
plementary information for a better understanding of the 
key role played by duality of the AA model. 



II. AUBRY- ANDRE MODEL 

The quantum dynamics of the AA model for the am- 
plitude ip n (t) is governed by the equation 

idtipn = -[^n+i + VVi-i] + Acos(27ran + ip)ij) n (1) 

where A is the disorder strength, and ip an arbitrary 
phase. This equation can be used to model non- 
interacting atoms subject to two periodic optical lattices 
with different wavelengths, that is, a bichromatic lattice 
[5, 6 . In this case, a is the ratio between their wave- 
lengths. In the tight-binding regime, Eq. is easily 
obtained starting from the Schrodinger equation and pro- 
jecting the continuous wave function ip(x) on the basis of 
the Wannier functions, w n (x), of the lowest band of the 
primary lattice [16], tp(x) = J2 n VVi^n^), where n is the 
lattice site index. The secondary lattice appears in the 
last term of Eq. u\l in the form of a modulation which 
mimics a disorder (quasi-disorder). Without any loss of 
generality, one can choose a < 1, because Eq. ([I]) is in- 
variant under a shift of a by an integer number. When 
a is a rational number one can write it as a — p/q and 
the solution of Eq. can be restricted to a region of 
size N = q, which coincides with spatial periodicity of 
the system. The case of irrational a can be obtained as 
a limit of a continued fraction approximation. 

Aubry- Andre showed [9] that this model has a duality 
under the following transformation 

h = 7V-V2 ^ n e in ^ l+ ^e-^ 1 (2) 

n 

which maps Eq. ([!]) into an equation for the new variable 
exactly of the same form as Eq. ([TJ but with disor- 
der strength 4/ A. This transformation corresponds to a 
projection on a basis of quasi-momentum eigestates with 
eigenvalues £ = 2ital + = (27r/N)pl + where £ is 
defined in the first Brillouin zone, £ G [— 7r,7r]. By using 
this duality one can show that the AA model undergoes 
a transition from an extended to a localized regime at 
A = 2. 



2 




35 

30 

25 

20 

15 

10 

5 



20 

15 

10 

5 





10 
Time 



15 



20 



FIG. 1: (Color online) Quasi-momentum distribution \fk{t)\ 2 
obtained from the DFT of the solution of Eq . The index £ 
is related to k by £ = 2nk/N. Here we use a = 0.2282... and 
(p = 0. The initial wave packet in real space is a Gaussian 
of width a = 10. Time is given in dimensionless units. Top 
panel: A = 1, only few modes are involved and a periodic 
oscillation of the central and side peaks is observed. The side 
peaks are at a distance ±27ra from the central peak. Bottom 
panel: A = 5, the evolution is affected by the coupling of 
many modes and the periodic oscillations are no more visible. 



In the AA model the space is discrete and one can con- 
veniently introduce the discrete Fourier transform (DFT) 



fk = N 



-1/2 



^2 A 



- kn 



(3) 



which corresponds to a projection on the basis of quasi- 
momentum eigenstates with eigenvalues £ = (2ir /N)k. 
The relation bewteen momentum and quasi-momentum 
distributions is \ip(k,t)\ 2 = N\f k (t)\ 2 \w(k)\ 2 where w(k) 
is the Fourier transform of the Wannier function centered 
on the lattice site n = 0. Note also that the expression 
(|3| is related to the dual mapping Q by an arbitrary 
shift of and a permutation [13]. By applying the trans- 
formation ([3| to the AA model one gets [16] 



idtfk = -2 cos 



/2tt/c 



fk- 



X 



-^h +p + e^f k - p ] (4) 



which is the equation describing the evolution in momen- 
tum space. Owing to the duality of the AA model under 
the transformation ([2| , and to the similarity of the latter 
with Eq. ([3|, the localization properties in momentum 
space are the same of the dual AA model, except for the 
fact that disorder couples modes differing by |A£| = 2na 
instead of neighboring ones. 



III. RESULTS AND DISCUSSIONS 

In the following we will consider the evolution of a 
wave packet, by solving Eq. ([I]) in real space and then 



computing the evolution in momentum space by means of 
the DFT [17 . As initial condition we choose a Gaussian 
wave packet, ipj(0) = Aexp{— j 2 /2cr 2 }, where a is the 
width and A is a normalization factor. This choice is 
convenient if one wants to simulate realistic experimental 
configurations; it also allows one to explore the behavior 
of sharp to broad wave packets in a continuous manner. 

According to the previous discussion, localization in 
momentum space occurs for A < 2, where the wave 
packet instead spreads in real space. In this regime one 
thus expects to see only one or few momentum com- 
ponents significantly populated. Conversely, for A > 2 
the regime is diffusive in momentum space and localized 
in real space, and one should see a momentum distri- 
bution with many modes coupled together during the 
evolution of the system. This is indeed confirmed by 
our numerical simulations, as shown in Fig. [I] (the value 
a = 1064.4/866.6 - 1 = 0.2282... has been chosen in or- 
der to model the bichromatic lattice of the experiment 
of Ref. [6]). A striking feature is that, for A < 2, the 
quasi-momentum components \fk\ 2 exhibit periodic oscil- 
lations, occurring among the central peak and two side 
peaks at distance ±27ra. The numerical result for the 
frequency of these oscillations is plotted in Fig. [2] as a 
function of A (red points for A < 2) [T8] . 

This oscillating behavior can be understood in terms 
of the following model. Let us assume that the width 
l/cr of the initial quasi-momentum distribution is small 
enough, so that only the k = mode can be considered 
populated at t = (/fc(0) = 5 ky o) and the time evolu- 
tion couples the mode at k — with two side modes at 
k = ±p only (i.e., £ = and ±27ra, respectively). This 
assumption is valid when l/cr < 2ixcx. In this way, the 
AA equation is mapped into an eigenvalue problem of a 
3x3 matrix, whose eigenvectors and eigenvalues can be 
written as g k j and Ej, respectively, with j = 1, 2, 3. The 




FIG. 2: (Color online) Oscillation frequency of the central 
peak of quasi-momentum (A < 2, red) and spatial (A > 2, 
blue) distributions, for a = 10, cp = 0. The full solution of 
the AA model (dots) is compared with the predictions of an 
analytic three-mode approximation (full lines) and a semi- 
analytic five-mode approximation (dashed line). 
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FIG. 3: (Color online) Spatial distribution \i/j n \ 2 obtained 
from the solution of Eq |T]) using a single site initial condition 
^ n — S n ,o. As in Fig. |lj we use a — 0.2282... and ip = 0. 
Time is given in dimensionless units. Top panel: A = 1, the 
initially localized wavepacket spreads ballistically and there 
are no visible periodic oscillations. Bottom panel: A = 5, 
the wave packet is localized and a periodic oscillation of the 
central peak and its nearest neighbors is clearly visible. 
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FIG. 4: (Color online) Visibility of the oscillations in real and 
quasi- momentum space a function A, for ip = 0. The numer- 
ical points are compared with the three-mode approximation 
(full lines). Bottom: only few modes have been initially pop- 
ulated by using a ^-function initial wave packet in real space 
and a Gaussian with a — 10 in momentum space. Top: many 
modes have been initially populated by inverting the initial 
conditions with respect to the bottom panel. 



initial condition is fk(0) = ^2j^\lj9k,j^ where the coef- 
ficients 7j are given by the standard rules of quantum 
mechanics. Under these assumptions one has 7^=3 = 0, 
and the time evolution takes the form 

l/fcW| 2 =(7i^M) 2 + (72^, 2 ) 2 

+ 7l72#/c,l#fc,2 cos [(^2 - E^t] . (5) 

This expression describes a time-periodic oscillation of 
the relative intensity of the central and side peaks, with 
frequency v(\ < 2) = \E2 — Ei\/2ir, given by 

v{\ < 2) = 7r- 1 v / [ 1 -cos(27ra)] 2 + A 2 /2. (6) 

It is worth stressing that, once a is fixed, this frequency 
depends only on the disorder strength A, but not on the 
phase cp. This three- mode approximation provides a rea- 
sonable description of the numerical results, as shown by 
the solid line for A < 2 in Fig. |2j 

The three-mode approximation becomes inaccurate 
when approaching A = 2, where more modes are cou- 
pled during the evolutions. In order to check this ef- 
fect, one can go one step fur her and consider a five- mode 
approximation in which the time evolution couples the 
modes k = 0, k = ±p and k = ±2p (i.e., £ = 0, ±27ra, 
and ±47ra). This is straightforward generalization of 
the three-mode approximation, except for the fact that 
the differential equations for the coefficients jj (t) do not 
yield simple analytical expressions and, moreover, the 
solutions contain several oscillation frequencies. The red 
dashed line in the A < 2 part of Fig. [2] is our numerical 
result for the dominant component of the frequency spec- 
trum, solution of the five-mode approximation, which 
mostly determines the time evolution of the central peak. 



As expected we get a better agreement with the full in- 
tegration of Eq ([I]) compared to the three-mode approx- 
imation, especially in the region close to the transition 
point A = 2. 

In the region A > 2 the few-mode approximation is 
expected to fail in the quasi-momentum space, where 
the wave packet is no more localized. Indeed, in this 
regime, we do not see any significant evidence of periodic 
behaviors in the quasi-momentum distribution (see the 
bottom panel of Fig. [T]). Conversely, owing to the dual- 
ity of the AA model, one expects periodic oscillations to 
take place in real space, where the wave packet is local- 
ized. This is confirmed by our numerical integration of 
Eq. 0, as shown in Fig. [3j In the top panel the wave 
packed spreads ballistically and one cannot detect any 
significant periodic behaviour; conversely in the bottom 
panel the wave packet is localized, while the central peak 
and its nearest neighbours oscillate periodically. By as- 
suming that the initial density distribution is localized 
in a single lattice site, n = 0, which is coupled with the 
nearest neighboring sites, n = ±1, we obtain a three- 
mode approximation analogous to the one used before in 
quasi-momentum space, but describing oscillations in the 
spatial distribution. The scenario in real space is more 
complicated because one generally observes oscillations 
with several frequency components, which also depend 
on the phase cp. However, in the special case ip = 0, one 
finds just a single frequency, given by 

u(\ > 2) = (27r)-VA 2 [l-cos(27ra)] 2 +8, (7) 

which is shown as the solid line for A > 2 in Fig. |2j In 
the same figure we also plot the frequency obtained from 
the full numerical integration of Eq. ([T]) (blue dots in the 
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A > 2 region); in this calculation we have used a Gaus- 
sian of width cr = 10 as initial shape of the wave packet, 
but we have also checked that the frequency v does not 
depend on a, except close to A = 2. The dashed line is 
the result of a straightforward semi-analytic extension to 
five modes, as in the A < 2 region. It is worth stressing 
that the condition for the validity of the few-mode ap- 
proximation for the oscillations in real space (ip n = S nj o 
or, equivalent ly, a < 1) is much more constraining than 
the one in momentum space (1/cr < 2ira) from the point 
of view of experimental realization. 

The amplitude of the oscillations in both real and mo- 
mentum space also changes with A, affecting its visibility. 
The latter can be calculated from the frequency spec- 
trum of the numerical solution of Eq. ([!]), as the ratio 
between the modulus of the Fourier component of fre- 
quency v(X) and the modulus of the component at zero 
frequency. In a consistent way, one can define the visibil- 
ity in the three-mode approximation; for the oscillations 
in momentum space for A < 2, the visibility can be writ- 
ten as (1/2)717200, 100, 2/[(7i#o,i) 2 + (720o,2) 2 ]- A similar 
definition can be given in real space for A > 2. In Fig. [4] 
we show the visibility of the oscillations as a function of A. 
The points are the numerical results, while the lines rep- 
resent the three-mode approximation. We have used two 
values for the width of the initial Gaussian wave packet, 
namely a = (i.e., a S- function) and a = 10. In the 
upper panel, the two values of a are used for A < 2 and 
A > 2, respectively. They correspond to a broad initial 
wave packet both in momentum space for A < 2 and real 
space for A > 2. In the bottom panel we use again the 
same values of a, but in the opposite regions, so to have 
a narrow initial wave packet in both spaces [19] . One 
can see that the visibility depends significantly on both 
a and A. Again, the three-mode approximation is quali- 
tatively correct, except near A = 2. We observe that the 
three-mode approximation gives a better agreement for a 
narrow initial distribution (lower panel), as in the oppo- 
site case of a broad distribution many modes are initially 
excited and this approximation is not expected to be ac- 
curate. Another interesting feature is the effect of the 
duality of the AA model. Indeed, in both panels, the re- 
sults in the region A < 2 almost coincide with those in the 
region A > 2 under the change of variable A — ^ 4/ A, pro- 
vided the initial distributions are broad (upper panel) or 
narrow (lower panel) in both momentum and real spaces; 
this duality also implies the continuity at A = 2. 



IV. SUMMARY 

We have shown that the time evolution of a wave 
packet in the Aubry- Andre model exhibits interesting pe- 
riodic behaviors both in momentum space, for A < 2, and 
in real space, for A > 2. The occurrence of oscillations in 
the momentum and density distributions can be used to 
test the applicability of the AA model to the description 
of the dynamics of ultracold gases in ID bichromatic lat- 




FIG. 5: (Color online) Phase- averaged intensity of the cen- 
tral peak in momentum space, \fo(t)\ 2 , as a function of time 
and disorder strength A for a wavepacket with a = 10. The 
intensity is given in arbitrary units. 

tices. We have numerically calculated the frequency and 
visibility of these oscillations and we have introduced a 
simple few-mode approximation to interpret their behav- 
ior. The use of momentum space has some important 
advantages. First, the frequency of the oscillations does 
not depend on the relative phase cp of the two lattices 
or, in other terms, on the initial position of the wave 
packet, which would be hardly controllable in typical ex- 
periments. Second, the condition for the applicability of 
a few-mode approximation is less restrictive than in real 
space, because the width of the initial wave packet in mo- 
mentum space can be easily made much smaller than the 
distance between the modes which are coupled in time 
evolution. 

Our analysis suggests that the oscillations of the cen- 
tral and side peaks in the momentum distribution can be 
efficiently used to probe the transition from diffusion to 
localization in the AA model. A possible strategy con- 
sists of measuring the intensity of the central peak as a 
function of time for different values of A, exploiting the 
fact that for A > 2 the oscillations are phase dependent, 
while for A < 2 they are not. Actually, in typical ex- 
periments with ultracold atoms, the phase (p varies at 
random at each realization, so that performing an av- 
erage over many realizations at fixed A is equivalent to 
an average over numerical simulations with different cp. 
Thus one expect that the oscillations vanish for A > 2 
(phase sensitive regime), but remain clearly visible for 
A < 2 (phase independent regime). This is shown in 
Fig. [5j where the average has been done over 50 different 
values of the phase cp for each value of A. Indeed the be- 
havior of \fo(t)\ 2 exhibits a transition at A = 2. From the 
same figure one can also extract the frequency ^(A < 2). 
By using the experimental parameters of Ref. [6 , with 
a = 1064.4/866.6 and A = 1, the oscillation period turns 
out to be of the order of 5 ms. 

These results are relevant for current experiments with 
ultracold atoms, where the momentum distribution can 
be detected with good resolution by performing time of 
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flight measurements. This study provides also a start- 
ing point for future investigations on the interplay be- 
tween interaction and localization. It is known that in 
certain regimes the presence of a repulsive interaction 
can destroy the disorder induced (Anderson) localization 
giving a sub-diffusive spreading of an initially localized 
wavepacket [lOj [11], but its experimental detection in 
real space is currently limited by the difficulty to reach 



the long time regime, and by the finite resolution in the 
observation of low density tails of the atomic distribu- 
tion [6 j. In this sense, the observation of oscillations in 
momentum space can be a more reliable tool. 
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